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1. Introduction 

Bilinear filtering or bilinear interpolation is used in computer graphics to com- 
pute intermediate values for a two-dimensional regular grid. One of the main ob- 
jectives is the smoothening of textures when they are enlarged or reduced in size. 
In mathematical terms, the interpolation technique is based on finding a function 
f(x, y) of the form f(x, y) = a + bx + cy + dxy, where a, b,c,d € R, that passes 
through prescribed data points. 

As textures reveal, in general, a non-smooth or even fractal characteristics, a 
description in terms of fractal geometric methods seems reasonable. To this end, 
the classical bilinear approximation method is replaced by a bilinear fractal in- 
terpolation procedure. The latter allows for additional parameters, such as the 
box dimension, that are related to the regularity and appearance of an underlying 
texture pattern. 

We introduce a class of fractal interpolants that are based on bilinear functions of 
the above form. We do this by considering a more general class of iterated function 
systems (IFSs) and by using a more general definition of attractor of an IFS. These 
more comprehensive concepts are primarily based on topological considerations. 
In this context, we extend and correct some known results from |12) concerning 
fractal interpolation functions that are fixed points of so-called Read-Bajraktarevic 
operators. Theorem [4] relates the fixed point in Theorem [3] to the attractor of an 
IFS and generalizes known results to the case where the IFS is not contractive. 

As a special example of the preceding theory we introduce bilinear fractal inter- 
polants and show that their graphs are the attractors of an underlying contractive 
bilinear IFS. Finally, we present an explicit formula for the box dimension of the 
graph of a bilinear fractal interpolant in the case where the data points are equally 
spaced. 
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2. General iterated function systems 

The terminology here for iterated function system, attractor, and contractive 
iterated function system is from |5J. Throughout this paper, (X, d) denotes a com- 
plete metric space with metric d = d\. 

Definition 1. Let M £ N. If f m : X — » X, m = 1,2, ...,M, are continuous 
mappings, then J- = (X; fx, /2, /m) *s called an iterated function system 
(IFS). 

By slight abuse of terminology we use the same symbol T for the IFS, the set of 
functions in the IFS, and for the following mappings. We define J- : 2 X —> 2 X by 

HB) := (J f{B) 

for all B £ 2 X , the set of subsets of X. Let H = H(X) be the set of nonempty 
compact subsets of X. Since T (H) C H, we can also treat J 7 as a mapping T : 
M -> H. When U C X is nonempty, we may write H(i7) = H(X) n 2 U . We denote 
by \F\ the number of distinct mappings in T . 

Let (in denote the Hausdorff metric on H, defined in terms of dx- A convenient 
definition (see for example p. 66]) is 

d m {B, C) := inf{r > : B C C + r,C d B + r}, 

for all B, C £ EL For 5 C X and r > 0, 5 + r denotes the set {y £ H : 3x £ S s.t. 
dx(x,y) < r}. 

We say that a metric space X is locally compact to mean that if C C X is compact 
and r is a positive real number then C + r is compact. The notation C + r means 
the closure of the union of balls of radius r, one centered on each point of C. 

The following information is foundational. 

Theorem 1. (i) The metric space (H, dm) is complete, 
(ii) If(K,d%) is compact then (H, dm) is compact. 
(Hi) Z/(X,(ix) is locally compact then (H, dj&) is locally compact. 

(iv) If X is locally compact, or if each f £ J- is uniformly continuous, then 
T : H — > H is continuous. 

(v) If f : X — >X is a contraction mapping for each f £ T , then T : H — > H is a 
contraction mapping. 

Proof, (i) This is well-known. A short proof can be found in p. 67, Theorem 
2.4.4]. 

(ii) This is well-known; see for example [9]. Here is a short proof. Let e > 
be given. Since X is compact we can find a finite set of points X £ C X such that 
X = U xe x e jB (x, e) where B(x,e) C X denote the open ball with center at x and 
radius e. Let H e := 2 Xe , a finite set of points in H. It is readily verified that 
HI = UceM c B {C, e) where now B(C,s) C II denotes the open ball with center at 
C £ HI and radius e, measured using the Haudorff metric. It follows that HI is 
totally bounded. It follows that HI is compact. 

(iii) Let C € HI. Consider the set C + r. It belongs to HI since X is locally 
compact. Let e > be given. Since C + r is a compact subset of X we can find a 
finite set of points C e C C + r such that C + rC U c£ c E B (c, e). Let C £ := 2 Ce , a 
finite set of points in HI. It is readily verified that C + rC U ce c c ^ (Cj £ ) where now 
B (c, e) C H denotes the open ball with center at c £ HI and radius e, measured 
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using the Haudorff metric. It follows that C + r is totally bounded. It follows that 
C + r is compact. 

(iv) Let B E H. We show that T : EI — > H is continuous at B. We restrict 
attention to the action of T on B + 1. If X is locally compact, it follows that B + 1 
is compact. It follows that each / € J 7 is uniformly continuous on 5 + 1. It follows 
that if X is locally compact, or if each / E J- is uniformly continuous, we can find 
S e > such that dx(f(x),f(y)) < e whenever dx(x,y) < 5 e , for all x, y G B + 1, 
for all f E T . Let C G H with dm(B, C) < S e and let f E J-. We can suppose that 
5 e < 1. 

Let 6' G /(-B). Then there is 6 G B such that /(&) = 6'. Since d w (B,C) < 6 e 
there is c G C such that d(b, c) < S £ . Since 5 e < 1 we have c € B + 1. It follows 
that d(f(b),f(c)) < e. It follows that f(B) c /(C 1 ) + e. By a similar argument 
f(C) C /(B) + e. Hence d M (f(B)J(C)) < e. 

(v) This is Hutchinson's theorem, [TOj p. 731], proved as follows. Verify that 
if A > is a uniform Lipschitz constant for all / G J 7 , namely d%(f (x) , f (y)) < 
Xdx(x,y) for all f E J 7 , for all E X, then A is also a Lipschitz constant for J 7 , 
namely rf H (^"(-B), -^(C)) < Ad H (B, C) for all B, C E H. If / : X ^X is a contraction 
mapping for each / G J 7 , then we can choose A < 1. It follows that T : H — > H is a 
contraction mapping. □ 

For B C X, let J rfc (B) denote the fc-fold composition of J 7 , i.e., the union of 
f il o / J2 o • • • o / i(s (B) over all finite words ■ ■ ■ ik of length k. Define ^(B) := B. 

Definition 2. j4 nonempty compact set A C X is said to 6e an attractor of the 
IFS T if 

(i) T{A) = A and 

(ii) there is an open set U C X such that A C U and lim/ c _ s . 00 J- k {B) = A, for 
all B E M(U), where the limit is with respect to the Hausdorff metric. 

The largest open set U such that (ii) is true is called the basin of attraction 
(for the attractor A of the IFS T). 

We will use the following observation jTTJ Proposition 3 (vii)], p. 68, Proposi- 
tion 2.4.7]. 

Lemma 1. Let {Bk}^ =1 be a sequence of nonempty compact sets such that B^+i C 
Bfc, for all k. Then nfc>iBj. = lim^oo B^ where convergence is with respect to the 
Haudorff metric. 

We use the notation 5* to denote the closure of a set 5*. 

Theorem 2. Let T be an IFS with attractor A and basin of attraction U. If J- : 
M(U)— > M.(U) is continuous then 

A = P| (J T k {B) for allBcU such that B G U(U). 

K>1 k>K 

The quantity on the right-hand side here is sometimes called the topological upper 
limit of the sequence {B fc (B)} fc _ . 

Proof. We carry out the proof under the assumption that B G M(U). (It then 
follows from [TT1 Proposition 3 (i)] that Theorem [2] is true for all B C U such that 
B G M(U).) 
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First we show that [J k>K F k (B) is a nonempty compact subset of U. It is 
nonempty because it contains F k (B) which is nonempty because B G H(Z7). Let 
{Oi : i £ 1} be an open cover of [J k>K T k (B). Then {Oi : i £ 1} is an open 
cover of A. Hence it contains a finite subcollection, say {O m : m = 1,2, --,M} 
such that A C U*f =1 O m . Since F k {B) converges in the Hausdorff metric to 
A it follows that, for sufficiently large K, we have T k (B) C U^ =1 O m for all 
fc > K x . It follows that Ufc^A'- 77 ^ 5 ) c u m=i°m for all if > ifi. It follows 
that, if K > Ki then [J k>K J~ k (B) is compact. If A" < Ki then we note that 
Uk>K J 7k (B) = {Ji< 1 >k>K J~ k (B) U {J k>Kl J~ k (B) which is a finite union of com- 
pact sets, and so is compact. We conclude that [J k>K F k {B) is a nonempty compact 
set. It is also a subset of U since each term in the union belongs to U, so the union 
does belongs to U, so the closure of the union belongs to U. 

Since f] k \J k > K+ iF k (B) C n k L) k > K T k (B) it follows that 



A = n u ^ fe ( s ) 



K>1 k>K 



is a well-defined nonempty compact subset of U . 
Now observe that 



T 



{A\ = T j ^lim |J F k {B) J (by Lemma [T]) 

y k>K J 

= lim J 7 II J rfc (B) (since T : H — > H is continuous) 

\fc>A / 

lim |J / I |J J*(B) (by the definition of F) 




is cts. 



lim U / I LJ ■7~ k (B) I (since is a finite union) 



K— >oo 

/GJ 7 \fe>if / fer 



= lim II II f(F k {B)) (by direct comparison, see below,) 
?r — ^ — ^ — 



fc>A /eJ 7 



lim II F k (B) = A 



A-s-oo 

fc>A 

(It is readily verified that / ^Jk>K^F k {B)) = \Jk>Kf{F k {B)) which implies the 
penultimate step above.) 

Since A is nonempty, compact and lies in U, the basin of attraction of A, we 
must have 

A = lim F k {A) = A. □ 



We will also need the following observation. 
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Lemma 2. Let X be locally compact. Let T — (X; /i, /a, /jv) be an IFS with 
attractor A and basin of attraction U . For any given e > there is an integer L 
such that for each x G A + e there is an integer I < L such that 

d H (A,P({x}) <e. 

Proof. For each x G A + e there is an integer l(x,e) so that du(A,J rl ^ c ' E '({x}) < 
s/2. 

Since X is locally compact it follows that lf l \ x ' e > : H — y H is continuous. Since 
jrl(x,e) . jj jj jg continuous there is an open neighborhood iV({x}) (in H) of {a;} 
such that d H {A,p( x - £ \Y) < e for all Y G JV({a;}). It follows, in particular, that 
there is an open neig hborhood N(x) (in X) of x such that d H (A, P<. x > e ) ({y}) < e 
for all y € N(x). Also since X is locally compact, there is a finite set of points 
{x%, X2, such that A + e c U^ =1 JV(xi). Choose L := max^ e). □ 

3. Fractal interpolants as fixed points of operators 

Let {(Xj, Yj) : j = 0, 1, iV} denote the cartesian coordinates of a finite set of 
points in the Euclidean plane, with 

Xq < X\ < ... < Xff. 

Let 7 denote the closed interval [Xo,Xn]. For n = 1,2, ...,N, let Z„ : L —> 
[X n -i,X n ] be a continuous bijection. Let L : I I be bounded and such that 

L(x) := C 1 ^) for x G (X„_i,X„) 

for n = 1,2, ...,iV. Let 5 : [Xo,Xjv] — > M be bounded and piecewise continuous, 
where the only possible discontinuities occur at the points in {Xi,X2, Xjv-i}- 
Let 

s := max{|S'(x)| : a; G [X , X N ]}. 

Denote by C — C(I) denote the set of continuous functions / : I —> K. It is 
well-known that (C, (ioo) is a complete metric space, where 

doo{f,g) = max{|/(i) - g(x)\ : x G /}. 

Let 

C7* := {/ G C : /(X ) - Y , f(X N ) = Y N }, 
C** :={feC: f{X 3 ) = for j = 0, 1, N}. 

Note that C* and C** are closed subspaces of C, with C** C C* C C. We say that 
each of the functions in C** interpolates the data {(Xj,Yj)}. 
Let b G C* and h G C7**. Define T : C C by 

T(s) :=fc + S-(soL-()oL). 

T is a form of Read-Bajraktarevic operator as defined in [15]. The following result 
is a corrected version of [12j Theorem 5.1, p. 136]. See also [10l Theorem 3, p. 
731]. 

Theorem 3. The mapping T : C — >• C obeys 

doo(Tg,Th) < sd^^g^h) 

for all g,h G C. In particular, if s < 1 i/ien T is a contraction and it possesses a 
unique fixed point f G C** . 
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Proof. To prove the inequality observe that 
dooCTfoJ.TCfc)) - max{\S{x){g(L(x)) - h(L(x)))\ : x e /} 

< smax{|( fl (|- 1 (a:)) - MC 1 ^)))! : * e n = 1, 2, N} 

= sdooig, h). 

The existence of a unique fixed point / € C (when s < 1) follows from the contrac- 
tion mapping theorem. Since f(C*) C C** and (C**, doo) is closed, hence complete, 
it follows that / € C**. □ 

Note that Tg = H + S ■ g o L where if = h — S -bo L. This tells us that a fractal 
interpolation function / is uniquely defined by three functions if, S, and L, of the 
special forms defined above. 

The fixed point / of T interpolates the data {(Xj,Yj) : j = 0, 1, 2, ...,N} and is 
an example of a fractal interpolation function lj. One way to evaluate / is to use 

/= lim T k (f ), 

where fo £ C* . The rate of convergence is governed by 

\\f-T k (fo)\\ 00 <s k \\f-f \\ 00 . 

4. The metric space (f x R, d q ) 

The following metric generalizes the "taxi-cab" metric. We will need it in the 
proof of Theorem [4] 

Proposition 1. Let a, /3 > and 3 : i ->• R. Let d, : (/ x I) x (/ x R) 4 [0, 00) 

6e defined by 

dqdx^y-s) , (2:2,2/2)) := a |a?i - x 2 \ + f3 \(y 1 - q(xi)) - (y 2 - q(x 2 ))\ , 

for all (2:1,2/1), (2:2,2/2) £ -fx R. Then d q is a metric on I x R. If q is continuous 
then (f x R, d q ) is a complete metric space. 

Proof. Clearly d q ((x 2 , 2/2) , (xi, 2/1)) = d 9 ((afi, 2/1) , (2:2, 2/2)) > 0. Suppose that 
d q {{x\, 2/1) , (2:2,2/2)) = 0. Then a|a;i -x 2 | + j8|(yi - 6(0:1)) - -6(a!2))| = 
which implies 2:1 = 2:2- Hence |(?/i — 6(2:1)) — (2/2 — 6(xi))| = which implies 

2/1 = J/2- 

Demonstration that if obeys the triangle inequality. Let (2^,2/i) € i X R, for 
i = 1,2, 3. Write = g(x^) for i — 1,2, 3. We have 

d<r((£i,2/i) , (2:2,2/2)) +d q {(x 2 ,y 2 ) , (2:3,2/3)) 

= a |xi - x 2 | + /3 I (2/1 - t?i) - (2/2 - <?2)| + a |x 2 - 2! 3 | + f) K2/2 - 92) - (2/3 - 93)1 
= aflaii — af 2 j + |aj 2 — a; 3 |) + P(\{yi - qi) - (2/2 + K2/2 - 92) - (2/3 - 93)1) 

> - x 3 |) -I- /3(|(2/i - qi) - (2/2 - 92)! + |(j/2 - 92) - (2/3 - 93)!) 

> a(|xi - 2; 3 |) + P{\{yi - qi) - (j/ 3 - g 3 )|) = dg((aci,j/i) , (2:3, 2/3))- 

To prove completeness in the case that q is continuous, let {(xk, 2/fc)}fe=i denote 
a Cauchy sequence with respect to the metric d q . Given e > we can find an 
integer N(e) so that 

a\x k -xi\+ (i\(y k - q(x k )) - (y t - q(xi))\ < e 

whenever k, I > N(s). It follows that {2:^} is a Cauchy sequence with respect to 
the Euclidean norm, and so it converges, with limit x* € I. Since q is continuous, 
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it now follows that {q(xk)} converges to some limit q* £ ML In turn, it follows that 
{Uk} converges to some y* £ E. Hence {(xk, yfc)}fcli converges to (x*,y*) £ I X K. 
It follows that (/ x M, d) is complete. □ 

5. Fractal interpolants as attractors of iterated function systems 

Here we characterize the graph of the fixed point / of T as an attractor of an 
IFS. Define w n : I xR ->J xlby 

w n (x, y) := (l n (x),h(l n (x)) + S{l n {x))(y - b(x))). 

Define an IFS by 

W := {I x R;wi,W2,...,w N ). 

Here we make use of the metric d q of Theorem [T] with q = /, the fixed point of 
T. Let B > and let 

X:={(x,y):x£l,\y-f(x)\<B}. 

It is readily verified that, when Theorem [3] holds, namely when s < l,yV(X) C X. 
The following theorem gives conditions under which (i) the IFS (X; wi,W2, ■■■,wn) 
is contractive with respect to df and (ii) W has a unique attractor. This result is 
a substantial generalization of |12[ Theorem 5.3, p. 140] which would require, in 
the present setting, that h is uniformly Lipschitz. Here, we avoid this restriction 
by using the metric d q with q = f . 

Theorem 4. Let s < 1 and let f £ C** be the fixed point ofT, as in Theorem^ Let 
l n : I — > I have Lipschitz constant Xi < 1, such that \l n {xi) — l n (x2)\ ^ \xi ~ 
for all X\,x 2 £ I, for all n. Let S : I — > [— s, s] have Lipschitz constant \$, so that 
\S(xi) — S(x2)\ < As \x\ — a?2 1 for all X\,X2 £ I- Then the IFS (X; u>\,W2, ...,wn) is 
contractive with respect to the metric df with a = 1 and < j3 < (1 — A;) /Ag-BAj. 
In particular, under these conditions, the IFS W has a unique attractor A = T(f), 
the graph of f , with basin of attraction I x R. 

Proof. Let (x 1: yi) ,(0:2,2/2) £ X. We have 

df(w n (xi,yi),w n (x2,V2)) - a \l n (xi) - l n (x 2 )\ 

= p\h(l n ( Xl )) + S{l n ( Xl )){ yi - b( Xl )) - f{l n ( Xl )) 

- h(l n (x 2 )) + S(l n (x 2 ))(y2 - b(x 2 )) - f(l n (x 2 ))\ 

= [3\ WlnfaVfa - f( Xl ))) - (S(l n (x2))(y 2 - f(x 2 ))) I 

< (3 \S(l n ( Xl ))\ -\( yi - f( Xl )) - (y 2 - f(x 2 ))\ 
+ \S(l n (x 1 ))~S(l n (x2))\-\(y2~f(x 2 ))\ 

< (3s\( yi - f( Xl )) - (y 2 - f(x 2 ))\ + (3XiX S B \x x - x 2 \ . 

Hence 

df(w n (xi,yt) ,w n (x2,V2)) 

< (aX t + flXsXiB) \ Xl -x 2 \ + ps\( yi - f( Xl )) - (j/2 - f(x 2 ))\ 
<(a + f3X s B) A ; \xi -x 2 \ + /3s\( yi - f( Xl )) - (y 2 - f(x 2 ))\ 
<c-(a\ Xl - x 2 \ +P\(yi- fix,)) - (y 2 ~ f(x 2 ))\) 

where c := max{s, A; + f3XsBXi/a} . Since A; < 1 we can choose a, /3 > so that 
c < 1. For example, we can choose a = 1 and < (3 < (1 — A;) /A5-BA;. 
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It follows that the IFS W := (X; w\, W2, • wn) is contractive, and hence it has 
a unique attractor. This attractor must be because a contractive IFS has a 
unique nonempty compact invariant set and it is readily verified that W(r(/)) = 
r(/). Since we can choose the constant B arbitrarily large, it now follows that W 
has a unique attractor, namely (But we have not provided a metric with 

respect to which W is contractive!) □ 

Note that 

T(T(g))=W(T(g)), for all g € C. 
When, for example, 5* is Lipschitz, s < 1, and the functions /„ are contractive, 
the graph of the fractal interpolant / can be approximated by the "chaos game" 
algorithm. (See [2] and [4] for new topological viewpoints of the "chaos game." ) 

6. Bilinear fractal interpolation 
We consider a specific example of the preceding theory. Let 

ln{x) ■■= X ^+{ X x~^X Q 1 ) (*-*0) 

and let 

S{x) :=s n {l-\x)) 
for x G [X n —i, X n ], n = 1, . . . , N, where 

s n (x) := s„_i + ( — - — y" 1 ) [x - X n _i) , 

\ A n — A„_l J 

with {sj : j = 0, 1, 2, TV} C (—1, 1). Then S is continuous and 

\S(x)\ <max{| Sn (CV))| :xe [X n _i, X n ], n e {1, 2, iV}} 
= max{|s.,| : j = 0, 1, TV} =: s < 1. 















-Vn — 





Let 



and let 

h(x) := y„_ a + I ) (.x - X„_x). 

\ A„ — A„_i / 

Theorem [3] implies that T has a unique fixed point /. Specifically, / is the unique 
solution of the set of functional equations 

f(l n (x)) - h(l n (x)) = [s n _i + (s n - s„_ 1 )x][/(a;) - 6(x)], x € /. 

We refer to / as a bilinear fractal interpolant. The reason for this name is that in 
this case the functions w n of the IFS W take the form 

w n (x,y) := (l n (x),a + bx + cy + dxy), 

where a, b, c, d are real constants. Functions of the form (x, y) M- a + bx + cy + dxy 
are called bilinear in the computer graphics literature. 
Specifically, 

(6.1) w n (x, y) = (k(aO,y n _i + %^ (x - X ) 



+ 



, + ( ■; B "- 1 )(»-^o) 



Ajy — A c 
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Figure 1 . A fractal interpolation function defined by three bilin- 
ear transformations. See text. 



In particular note that 

w n (X N , y) = {X n ,Y n + s n (y - Y N )) and w n+1 (X ,y) = (X n , Y n + s n (y - Y )). 

It follows that the images of any (possibly degenerate) parallelogram with vertices 
at (X ,Y Q ± H) and (X N ,Y N ± H ) , for H e K under the IFS fit together neatly, 
as illustrated in Figure [T] 



7. BOX DIMENSION OF BILINEAR INTERPOLANTS 

For present purposes we define the box-counting or box dimension of a bounded 
set M C R" to be 

(7.1) dua B M:= lim lo f^ M \ 

e^0+ loge^ 1 

where Af e (M) is the minumum number of square boxes, with sides parallel to the 
axes, whose union contains M. By the statement "dims M — D" we mean that the 
limit in equation (7.1) exists and equals D. 

Remark 1. It suffices to compute the box dimension of M using only discrete 
values of e, such as e := N~ m , for any KJVeN and m € N. 

In the case where M is the graph r(/) of a function /, knowledge of the box 
dimension of T(/) provides information about the smoothness of / since dims T(/) 
is related to Holder exponents associated with /. (See, for example, [131 Section 
12.5].) 

The next result gives an explicit formula for the box dimension of the graph of 
a bilinear fractal interpolant. The proof is based on arguments first applied in [7]. 
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Theorem 5. Let VV denote the bilinear IFS defined above, and let T(f) denote its 



attractor. Let a n — 1/N for n — 1,2, 
not a straight line segment then 



, N, and let ^2 



N 



-l+s„ 



> i. ifr(f) ™ 



log 



dim B r(/) = 1 



N 

\ - s r . 



n=l 



logiV 



otherwise dims L(/) = 1. 



Proof. Note that in the computation of the box dimension of L(/) it suffices to 
consider covers of T(f) whose elements are squares of side N~ r , r € No. Denote by 
Co(r) a cover of consisting of a finite number of squares of side N~ r , r £ No- 
Now consider a specific cover C(r) of T(f) of the form 



(7.2) C(r) := 



k - 1 



k 



a, a 



1 



r e N ; k = l,...,N r ; a e 



By the compactness of r(/), there exists a minimal cover Cg(r) of and also a 



minimal cover C*(r) of T(f) of the form ( 7.2 ). Denote by Afo(r), respectively, Af(r) 



the cardinality of these minimal covers. Since covers of the form (7.2) are more 



restrictive, we have Afo(r) < Af(r). On the other hand, every (N~ r x 7V _r )-square 
in Cg(r) can be covered by at most two (N~ r x iV~ r )-squares from a cover of the 
form (7.2). Thus, Af(r) < 2Afo(r). Hence, when computing the box dimension of 
r(/) it suffices to consider covers of the form (7.2). 

By an affine transformation, one can always achieve that yo = vn = 0. Since the 
box dimension is invariant under an affine transformation, we will restrict ourselves 
to this particular case. 

To this end, let r € No be fixed. Let C(r) be a minimal cover of of cardinality 
Af(r) consisting of squares of side N~ r whose interiors are disjoint. Let C(r,k) 
be the collection of all squares in C(r) that lie between x = 
k = 1, . . . , N r . Denote by Af(r, k) the cardinality of C(r, k), and let 

TZ(r,k):= |J C % . 

As C(r) is a cover of of minimal cardinality, every square in C(r) must meet 
L(/), and since / is continuous on /, the set lZ(r, k) must be a rectangle of width 
N- r and height h(r) := N-' r Af(r,k). Note that Af{r) = J2k=i Jf(r,k). 

Now apply the mappings Wi, i = 1,...,N, defined in ( |6.1[ ) to the rectangle 
lZ(r,k). The image of lZ(r,k) under Wi is a trapezoid contained in the strip 
l(k,i) - 1 l(k,i)~ 



and x = 



_/yr+i ' jyr+i 



x E, with l(k,i) :=k + (i- l)N r . Observe that 



N N r 



i=l k=l 



N 



The fixed point equation for L(/), namely, L(/) = Wi(T(f)), implies that 

i=i 

r(/)C \J Wi i\jTZ(r,k) . 



i=l \fc=l 
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Now, each trapezoid Wi(lZ(r,k)) is contained in a rectangle of width N and 
height 

— +fj, ik h{r), 

where \x ik := max{si_i + (sj-Sj_i)^i, Si-i + fa — Si-i)^?}. Note that s l _i + (s i - 

s»-i)^f- > and Si_i + (sj - Sj-i)^ > 0. We write /Lt ife = Si_i + (s t - Sj_i)^P, 
with e(fc) being either fc — 1 or fc. Therefore, 

W(r + 1, /(fc, *)) < + Mb ft(r)) N r+1 + 2 = Nfi ik M(r, fc) + {N\ ai \ + 2). 
Summation over k = 1, . . . , N r and i = 1 . . . , N produces 



Af(r + 1) < N 



N 



sn - so 
N r 



]Te(fc)AA(r,fc) 



k=l 



ci AT 



r+l 



where we set c\ := YliLx 



To estimate the sum over fc, we employ Corollary 3.3 in [5] with p k := 



and m :— 1. This yields 

N r 



_ jV(r,k) 
jV(r) 



N T 



s(k)Af(r, fc) < E kAf ( r > fe ) 



< 



k=l 



k=l 



.fc=i 



< 



A/"(r) 



k=l 



7v" + i)= A/ M(^ 



Therefore, 



A/"(r + 1) < N 7 + 



|sat - s | 



2AT 

= (iV 7 )[l + / 

where we set 7 := J2"=i s "~ 1 2 i ~ s ™ and /? := 
r finally yields 

■r-1 

3=0 



AT(r) + Cl N r+1 
(AT 7 )[1 + PN~ r ]j\f{r) + Cl N r+ \ 



- |s " 27 So1 ■ Note that 7 > 0. Induction on 



+ ciiV r 1 + 7 [1 + /37V- r ] + ■ • • + 7 r - x JJ [1 + pN- {r ~ 3) ] 



r I r — 1 T 

=(^7) r n[ i+ ^~'w) +ci ^ r i+ e^ n [i+z 3 ^] 

;=0 V 3=1 !=r- 1-fc 

Next, we show that the two products in the above expression converge absolutely as 
r — > 00. We consider the first product, the argument for the second one is similar. 
Note that {IlLoI 1 + PN^ 1 ] : r E N} converges absolutely iff {E[=o P N ~ l '■ r e N i 
converges absolutely. The latter sum, however, converges to Nf3/(N — 1). Denote 
the value of UZol 1 + P N ~ l ] h y P (/ 3 ) and observe that P(/3) < exp(7V/3/(iV - 1)). 
Depending on the value of 7, two cases need to be considered. 
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Case I: 7 < 1. This implies that Af(r) < N r [P(p)Af(l) + c x (l + (r - 1)P(/?))]. 
Hence, 

< lim log(^TOW)+ C i(l + (r-lW))]) = L 

r-s-oo logiV r 

Case II: 7 > 1. Here, we obtain 



r-l 

vJ'— r 



Thus, 



jV(r) < ( 7 A r ) r P(/3) + ci ( 7 iV) r \^~ r + P(P) J2 ^ 

< ( 7 AT) r [P(/3)AA(1) + d max{l, P(/3)}(1 - 7)- 1 ] =: c 2 ( 7 iV) r . 



^ r logC 2 {jN) r log 7 

dims < lim — — — = 1 



r-s-oo log N r log N 

To find a non-trivial lower bound for dims note that since / is a continuous 
function, dinisr(/) > 1. If T(f) is a line segment, i.e., if the set of data J :— 
{(Xj^Yj : j = 0, 1, . . . , N)} is collinear, then T(f) = [0, lj. Hence, dim B T(f) = 1. 

Suppose then that 7 > 1 and that is not a line segment, i.e., J is not 

collinear. Since each C; € C(r, fc) meets T(/), the image of Cj under the maps tUj, 
i = l,...,iV, must also meet r(/). Hence, C(r + l,l(i,k)) must at least cover a 
rectangle of width 7V~ (r + 1 ) and height (|^fc|[7V(r, fc) - 2] - |ai|) N~ r , where v ik := 
min{si_i + (si — s^i)-^ 1 , s,_i + (sj — Sj-i)-^}. This then implies that 

M{r + 1, i)) > N(\u ik \ [7V(r, fc) - 2] - H) - 2. 



Summation over fc and i, and using the rearrangement inequality 8, Equation 
10.2.1] together with the fact that if the sequence {Af(r,k) : k = l,...,N r } is 
arranged in decreasing order, 

ATr _ 1 " jy-r _ i 

fe=l k=l 

gives 

Af(r + 1) > (7JV) \l + t3N-( r+1 A Af(r)-c 3 N r+1 > N) Af(r) - c 3 N r+1 , 

N 

with C3 := 2 + -/V'y^ (|aj| + 2min{sj_i, Sj}). Induction over r yields 

i=i 

Af(r + 1) > ( 7 iV) m AA(r-m + l)-c 3 iV , ' +1 (l+7+---+7' n - 1 ) 

> ( 7 iV) r - m+1 AA(m) —. r 7V r+1 

1 — 7 

= ( 7 iV) r - m+1 (V(m) 



-1 



1-7 

for all m € N with 1 < m < r. 

To complete the proof of the theorem, the following lemma is needed. 
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Lemma 3. Suppose that f is a bilinear fractal in terp olant. Denote by Af(r) the 
cardinality of a minimal cover ofT(f) of the from (7.2). If 7 := X)i=i St ~ 1 2 ^ Sl > 1 
and ifT(f) is not a line segement then 

N r 

hm —T-- = 0. 

r-s-oo TV (r) 

Proof. The assumption that T(/) is not a line segement implies the existence of at 
least one index io € {1, . . . , N — 1} so that 

S:=\y i0 \ >0. 

Since / is continuous on [0,1], we have that Af(r) > SN r , for any finite code a 
of length r € N. Note that I is mapped to the line segments (a?i— 1, t/j-i), (xi,yi), 
implying that 

N 

Af(r) > + («i - s *-i) «o] * AT, 

i=l 

where we set xq '■= io/N. Inductively, this yields 

JV k 

■A/» > n [s ie -i + (s ie - Sit-i) L ie+1 ... ik (x )] 6N r , 

,i k =it=i 

with Li 1 ...i m := Ljj o • • • o Li m and -^(xq) := xq. We note that 



L il ... im (x) = ^ rn lx + Y / N n - 1 (i n -l) 



\ "=1 / 

for all m € N and x € [0,1]. The rearrangement inequality [El Equation 10.2.1] 
implies that for all to <E N 

m m 

^iV^n-l^EAr™-^*-!), 

n— 1 n— 1 

where the i* n £ {1, . . . , -/V} are taken in decreasing arrangement. Evaluating the 
latter sum and simplifying, one obtains the existence of an too € N such that 

Lii-i m (?o) > Ij Vm > to . 

Hence, 

■Af(r) > C7 r <5iV r , Vr > to > to . 
and for some constant c > 0. As, by assumption 7 > 1, the statement follows. □ 

Now, by Lemma [3j one can choose r and to large enough so that 

N{m) - ^— > 0, 

1-7 

and, therefore, Af(r) > €4(7 N) r , for a constant C4 > 0. Hence, dimsr(/) > 
log TV 

Remark 2. TTie proof of Theorem [3] shows in particular that for a given word 
a = a 1 . . . o r of finite length \a\, there exist constants < c < c such that 



c(7iV) |CT| <Af(\a\) < 5(7 iV) 



M 
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Moreover, if w ai ... lTr (T(f)) denotes the image ofT(f) under the maps w (Tl ... (Tr := 
w (Tl o ■ ■ • o w rTr over the subinterval L r7l ... rTr (I), then there also exist constants < 
c* < c* such that 

(7.3) c* X ai ■ ■ ■ A CTr iVH < J\f r7l ... rTr (\a\) < c* A CTl • • • A CT „ , 

where M ai ---a T d enote s the minimum number of N~\ a \ x iV - '""' -squares from 



a cover of the form (7.2) needed to cover w ai ...,j r (T(f)) and \ := s ' t Sv , i = 
1 V. 
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